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Abstract 

The quaternion Fourier transform (qFT) is an important tool in multi-dimensional data analysis, 
in particular for the study of color images. An important problem when applying the qFT is the 
mismatch between the spatial and frequency domains: the convolution of two qnaternion signals does 
not map to the pointwise product of their qFT images. The recently defined ‘Mnstard’ convolution 
behaves nicely in the frequency domain, but complicates the corresponding spatial domain analysis. 

The present paper analyses in detail the correspondence between classical convolution and the 
new Mnstard convolution. In particnlar, an expression is derived that allows one to write classical 
convolution as a finite linear combination of suitable Mustard convolutions. This result is expected 
to play a major role in the further development of quaternion image processing, as it yields a formula 
for the qFT spectrum of the classical convolution. 
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1 Introduction 

Since more than a decade ago, the use of quaternions has been investigated for 3D data analysis. As a 
guiding example we use color images. Indeed, one can represent, say, a color image as a pure quaternion 
via the identification 

R{x, y)i + G(x, y)j + B{x, y)k 

where x and y denote the pixel coordinates, and R{x,y), G{x,y) and B{x,y) are the red, green and 
blue color channels. When pixels are represented in such a way, they can subsequently be processed in a 
non-marginal way. This means that all color components of the image are treated together, as opposed 
to processing each of the three components independently. The main mathematical tool applied in this 
context is the quaternion Fourier transform (qFT) [TTJ [5T1 [5J [inj [TB]. This transform generalizes the 
classical Fourier transform by replacing the imaginary unit by two pure quaternions /i and v that square 
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to -1. These two quaternions can be chosen independently, but are typically related to relevant directions 
in the color space, such as the gray line or a direction perpendicular to it. 

So far, the quaternion approach has led to many successful applications. We mention edge detection 
and image filtering uni [9], flow visualization [lOj . watermarking [I], pattern recognition HU EH] , gener¬ 
alized wavelets [ZIEIII], noise removal from video images [16] and many more. In particular for saliency 
detection, which is the computerised detection of image features likely to attract human attention, use of 
a qFT method based on its phase spectrum vastly outperforms other methods in terms of efficiency and 
robustness M- 

The success of classical complex signal analysis is almost entirely due to the interplay between the 
Fourier transform (FT) and the classical convolution product and its application to the study of linear 
time-invariant (LTI) systems. Design of filters is very flexible and can be done both in the time and 
frequency domains. A fast computer implementation is obtained in the frequency domain, as the FT 
maps the computationally involved convolution of two signals to the much simpler pointwise product of 
their respective Fourier representations. 

The main issue in applying the qFT, as already indicated in |19j . is the mismatch between spatial 
and frequency domain: the qFT of the convolution of two quaternion signals does not equal the product 
of their respective qFTs. This means that, although e.g. the qFT saliency detection method works 
extremely well, it is not understood why it works so well. In turn, this means that other applications 
are not easily derived. The situation is even worse for the quaternion color edge filter |20j . which, by 
exploiting the geometric properties of the quaternions, generalizes the Prewitt edge detector and detects 
sharp changes of color in images. This filter is still implemented using the classical convolution product, 
thus leading to a slow algorithm, as the qFT cannot yet be applied in a meaningful way. 

This brings us to the aim of the present paper, namely to correctly connect spatial and frequency 
domain for the qFT, i.e. we want to match quaternion filter design in the spatial domain (using classical 
convolution) with design in the quaternion frequency domain (using pointwise products). This will allow 
us, in a series of follow-up papers, to construct a complete signal processing theory for quaternion LTI 
systems, which would be a huge breakthrough. 

The key idea to solve this problem is to find a suitable decomposition of the two signals to be convolved, 
so that the action of the qFT simplifies to a sum of products. This will be tackled using a mathematical 
tool we call Mustard convolution, introduced in [B] as follows: 

/ g := 

with / and g two quaternion signals and the qFT. We will establish a precise connection between 
classical convolution and Mustard convolution by expressing the standard convolution * as a finite linear 
combination of Mustard convolutions in the following way: 

f*g = Y.{aJ^^)*^,Ah9^n, eH. (I) 

i 

where denote suitable reflections on the arguments of the signals / and g (see Notation 12.41) . 

Subsequently applying the qFT to this formula then yields the desired connection between spatial domain 
and frequency domain, see e.g. Theorem 15. Il for the case of perpendicular roots. 

The expression m will be derived in three subsequent steps, depending on the relation between 
and V, the quaternion roots of -1 used to construct the qFT. We will distinguish 

• the case where g, = see Theorem 13.II 

• the case where g is perpendicular to see Theorem [331 

• the general case, see Theorem 14.II 

Whereas the case g = v \s relatively simple to deal with, the proof of the general case is hard. This 
lies in the fact that the quaternions do not form a commutative algebra and extra care has to be taken 
in order to swap exponentials, the roots g and v, and the quaternion signals / and g in order to arrive 
at the desired result. As the resulting theorems are complicated for arbitrary choices of roots, we have 
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implemented and thoroughly verified all our results in Theorems 12.5113.1113.31 and 14.11 using the symbolic 
algebra software Maple. 

The paper is organised as follows. In Section [5] we give basic definitions for quaternions and the qFT 
and derive several important identities we will need later on. In Section [3] we establish our result for 
equal and perpendicular roots. Next, in Section |4] we treat the general case of a qFT defined using two 
arbitrary roots of -1. Finally, in Section[5]we discuss some consequences of our results. 

2 Quaternions and the quaternion Fourier transform 

2.1 Quaternions 

The algebra of the quaternions HI is defined as follows: 

Definition 2.1. The quaternion algebra HI is the algebra over K generated by the basis elements i,j,k, 
under the relations 

= ijk = —1. 

These relations immediately imply: 

ij = —ji = k, jk = —kj = i, ki = —ik = j. 

A general quaternion q is characterized by four real numbers: 

q = a + bi + cj + dk. 

Sometimes it is useful to split a quaternion into a scalar and a vector part. The scalar part is the real part 
of the quaternion, S{q) = a. The vector part of the quaternion consists of the purely quaternionic part 
V (q) = bi + cj + dk. The quaternion product of qi = ai + bii + cij + dik and 52 = 02 + + C 2 J + c? 2 fe 

is given by 

9192 = 0102 — 6162 — C1C2 — did 2 + (0162 + bia 2 + cid 2 — diC 2 )i + (aiC2 — 61^2 + 0102 + dib 2 )j 
+ (01^2 + biC 2 — C162 + dia 2 )k. 

In the rest of the paper, we will always consider quaternion functions. These are functions / from 
to H. They can be expanded as 

f = fo + fii + / 2 J + /afc, 

where we will assume fi G L^(R^) so that Fourier transform and convolution are always well-defined. 

2.2 The quaternion Fourier transform 

In this paper we consider the left quaternion Fourier transform which is defined as follows: 

Definition 2.2. Denote by G I 2 = {a G Hja^ = —1} two geometric square roots of minus one, i.e. 
fjf = 12 ^ = —1. The left qFT of a function / : R.^ — >■ HI is defined as: 

= ^ f e-™e-™f{x)dx (2) 

with X = {xi,X 2 ) to simplify notations. 

The inverse transform is given by 

{f){y) = ^ [ e™e^^^^y^f{x)dx. 

271 Jr2 

There exist a plethora of different Fourier transforms in quaternion analysis, and the above definition 
is obviously not unique. Due to the non-commutative character of the quaternions the transform can 
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also be defined by sandwiching the function / with two exponentials, or by putting the exponentials on 
the right side. These are not the only three options, for an overview of different qFTs and their basic 
properties the reader can for example consult mm- In the present paper, we will always work with the 
left qFT, but the computations can be adapted to other definitions without additional difficulties. 

Based on an idea of Mustard in the context of the fractional Fourier transform [18], we define a 
generalized convolution for the qFT following | 6 ]. 

Definition 2.3. For the qFT , the Mustard convolution is given by 

if 9)ix) ■■= 2^{F^n-\^^''^if):F^'''i9))ix)- 


It clearly satisfies 

^^'''if *f.,.9) = 27TFf^’'^if)F>^’''i9). 

We introduce the following notation for reflecting the arguments of functions. 

Notation 2.4. For a function / : HU and a multi-index <p = {(pi, P 2 ) with (pi, (p 2 € {0,1} we put 

:= /((-l)^^xi, {-l)^^X2). 


Suppose furthermore that we are dealing with two square roots of 
anticommutator: 


a = {pL, v} = fiv vfi G K.. 


— 1, /r and n. We introduce their 


This anticommutator is always in R. Indeed, when pi and v are square roots of —1 they are pure 
quaternions (i.e. their real part is zero) and simple calculations reveal that they satisfy the following 
conditions: 

pL = hii T Cij -\- dik, 62 T T d-^ — 1 

V = b2i + C2j + d2k, 62 + ci + <^2 = 1 

a = {pL, v} = pLV VpL = —2(6162 + C1C2 + ^1^2), 

where i,j,k are the standard quaternion units. In particular, when pi and v are perpendicular, then 
a = 0 . 

We may now formulate the following theorem, see | 6 ], which expresses the Mustard convolution as a 
finite linear combination of classical convolutions: 


{f*9)ix)= f{y)g{x-y)dy. 

Theorem 2.5. The Mustard convolution can be expressed as, with {pi, = a G R., 


/ 9 = 


1 

4 1 
+ 


f *9 + f * + f *9 


( 1 . 0 ) 


di.i) 


—i/f * vg vf * — vf * lyg''^’^^ i/f * i/g 


(1.0) 


di.i) 




1 

+ - 


* pLvg + * pLiyg^^’^^ + * pLvg^^'^^ 

^ * yvg - f * - f * p.ug^'^'^l + / * p.vg^^'^1 

+ (fyf * yg + vf * -vf * -vf * 




1 ) 


* F9- * F9 


(1.0) 


* yg 


( 1 . 1 ) 
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Remark 2.6. The formula in Theorem \2.5\ simplifies dramatically for special choices of the roots p and 
v. When p = 2iv only four terms remain. When p and v are perpendicular (a = {p,v} = 0) 16 terms 
remain. 

Remark 2.7. It is possible to derive similar theorems for the case of right qFT, resp. the two-sided qFT. 
For the right qFT this leads to an equal number of terms. However, for the two-sided qFT the result 
simplifies and only 16 terms are necessary for the most general case of two roots p and v. 

The main aim of the present paper is now to obtain the ‘inverse’ result, namely to express the standard 
convolution * as a finite linear combination of Mustard convolutions in the following way: 

f *9 = a^, 6* G H. 

i 

As we will see, the proof of that statement is quite a bit more involved. 

2.3 Basic formulas 

In this section we summarize some of the basic formulas we will use to derive our results. Denote again 
by I 2 the set {a G ]HI|a^ = —1} of geometric square roots of minus one. By p, v we denote elements of I 2 
and hy q,r elements of El . 

2.3.1 Commuting and anticommuting of quaternions 

Following [ISIII], we define the commuting (j = 0) and the anticommuting (j = 1) part of a quaternion 
q with respect to p, a square root of —1 as follows: 

9wiu) = ^( 9 - (-1)V9M) 

q = qco ifj.) + qc^^L) = g- + ?+ 
q-p = pq- 
q+p = -pq+, 

where we introduced the notations 9cO(/i) = q- and qc^(p) = q+. These shorter notations are sometimes 
used when the quaternion we wish to commute is clear from the context. 

Now suppose we are dealing with two square roots of —1, p and v. Their anticommutator satisfies 
a = {p, v} = pv + G K. We can use this fact to simplify many calculations. The (anti)commuting 
parts of p with respect to v are for example: 

AicO(^) = p_ = ]^{p- vpv) = 

1 , . a 

Fc^{u) = T+ = i^KF + vpv) = p-\--v 

These commuting parts have interesting properties: 

P^=pI + P^_ = -1 (3) 

p+P--\- P-P+ = 0, ( 4 ) 

which can be derived by evaluating —v in two different ways. Indeed, we have 

—V = p^v 

= {t+F 

= {p\ + P+P- + P-P-\- + pt.)v 
= v{p\ - p+P- - p-p-^- + P-), 
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but also 


— V = I'll 

= + ^J^.) 

which leads to © and (IH) by equating both computations. 


2.3.2 Commuting and anticommuting of quaternions with exponentials 

In the study of the qFT we are mainly interested in commuting quaternions with exponentials such as 
= cos0 + ^sin 6 *, with = —1. We can derive the following commutation relations for such functions: 


9 = 9 -+?+= qcO{p,) + 

= e=F%±, 9eR. 

We will also frequently use its commuting and anticommuting parts: 

1 


( 5 ) 




= 2(6 - (-l)V(cos( 6 »)-z/sin( 6 l))/r) 
= + (-ly + a(-l)Vsin ( 0 )) 


( 6 ) 


where we recall that {^, v} = a. For the special case a = 0, which we will encounter later on, the formulas 
reduce to: 

1 




J COS0, if j = 0, 
I —i^sin0, if j = 1. 


( 7 ) 


In this case the commuting part is the cosine, while the anticommuting part is the part with the sine. 
We can apply our result with the purpose of swapping an element of I 2 and an exponential: 


He 


-vB 


-uB 


= (Ai+ +M-)e 

= e"V+ + e-"V- 


( 8 ) 


2 ' 2 

Later on we will also need the swapping of two elements of I 2 and their exponentials: 

= [e>^%+ + [e"V+ + e’^V-] 

= e''^[{v+)cP(y)^i+ + 

+ + {v+)c0{v)9-\ + + (i/_)cO(,,)/r_]. 

Notice that in this case the plus/minus-notation (/i+)+ would have been ambiguous. Easy calculations 
yield the following: 

(z/+)cO(i,) =v+-n- 


(*^+)ci(i') — 2^^+ 

/ \ _ a 


( 9 ) 
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Making use of dH]), and the property ([3]) and (|1]), we can eventually simplify the expressions to obtain: 




( 10 ) 


2.3.3 Properties of the quaternion Fourier transform 

We can translate the calculations on the exponentials to properties of the qFT as dehned in Definition 
12.21 A first important formula concerns the change of signs of the roots: 






( 11 ) 


which is obtained by performing a substitution in the integral (x —>■ x‘^) and using Notation 12.41 

In the case where ^ and ly are anticommuting it is also easy to obtain that passing /r and/or ly through 
the transform changes the sign of the corresponding exponential and thus of the arguments: 


27 !" Jr2 


( 12 ) 


In the case where {/x, v} = a, we can derive a similar property making use of equations (|S]) and (1101) . 
This yields after some computations: 


iy:F>^^''[g){u) = .F'^-'^(i/g(i’0))(u) + 

(13) 

as well as 

^ ^ (14) 


3 Important special cases 

Before dealing with the general case, we consider two important special cases. 


3.1 The case of equal roots 

When we have just one root g, = v the transform m reduces to: 

277 Jr2 


7 


Now we can compute the interaction with classical convolution by applying the definitions and subse¬ 
quently taking the integrals to the outside: 


* g{x)){u) = ^ * g){x)dx 

Jr2 

= 1 /■ f _ y)g(jj)dydx 

yR2 yR2 

= ^ f [ e-^^^^-^+^^-^^f{x-y)giy)dydx 

yR2 7 r 2 
27r Jr2 Jr2 
27'" 7r2 JrS 

In the last step we have made the substitution x ^ z + y. Because there is just one root, all exponential 
functions commute, and the exponential property holds. This will no longer be the case later on. We can 
continue by separating the parts in z and the parts in y. 


F^{f * g{x))(u) 
^77 JR2 JR2 

Z ZTT J^2 Jr2 

2 27r 7 r2 Jr2 V 


97r 

= -[:F>^{f){u)T^{g){u)-T^^{^^f){u)T>^{gg){u)+T>^{f){u)F>^{g<^^^^^){u)+T^{gf){u)F^{gg^^^^^^^ 

In a final step we apply the inverse Fourier transform to both sides. By using Definition 12.31 of the 
Mustard convolution and simplifying the notation f *fj, g ■= f g, we rewrite the result as 

/ * 5 = ^[/ y - m/ */x W + / *M + m/ *fj. 

This gives us our first result, which already showcases some of the important steps that will be made 
later on. We summarize it as follows. 


Theorem 3.1. Let f and g be quaternion functions on R^, and the left qFT with one root g. The 
classical convolution can then be expressed as a sum of Mustard convolutions by the following formula: 

f * 9 = ^[f *u 9 - Ff *u F9 + f *u 9^^’^^ + Tf *u T9^^'^^]- 

Remark 3.2. Note that one can consider in a very similar way a transform where the exponential factors 
are permuted in a different way, or with different signs in the exponential. 


3.2 The case of perpendicular roots 

Note that in this case a = 0 and the roots fi and v anticommute. We begin with the same steps in this 
more general case: 




*g){x)dx 


I 


-^XiUi —I/X 2 U 2 


fix - y)giy)dydx 


'R2 JR2 


e-^^^'^^e-’^--'^-fix-y)giy)dydx 


-[ [ 

27r 7 r2 7r2 


„-/j(zi+yi)«i p-i^(z2+y2)u2 


fiz)g{y)dzdy, 


where we again made use of the substitution x ^ z + y. It is now our purpose to rewrite the integrand 
such that all exponentials with 0 as integration variable are paired with the function / and those with 
y as integration variable with g. To do so, we will use the commutation relations for quaternions and 
exponentials as established in Section [231 

First, by the use of formula ([5|) we obtain for the two factors we want to interchange: 


-VZ2U2 _ „-izz2U2 -fj,yiui 


-J^Z2U2 ^{XyiUi 


(15) 


In the case of g and v anticommuting we can subsequently make use of the formula 0. We then continue 
in the following way: 


Z'X 7r 2 yR2 


g-M(zi+j/i)«i e-''^^^+y^^^^f{z)g{y)dzdy 


27T 


^-1>Z2U2 I pJ^Z2U2] 




+e 


'R2 J^2 

\^-l^Z2U2 _ ^1^Z2U2'\ 


^f.y^m^-i.y 2 U 2 f(^z)giy)dzdy. 


Now we still have to pull the function / through the exponentials. We obtain in a similar way as before: 






^-1^Z2U2 I ^,^"'22^21 


'R2 JR2 


[((/(^))c0(i'))c0(/z)6 


— liyiui —vy2U2 


iif i^))c°{v))c^{fi)^‘ 


izyiui p — vy2U2 


“b((/(^))ci(i'))c0(^)6 

\q-VZ2U2 _ g!^Z2«2l 


fiyiui vy2U2 


+ 9 iy)dzdy 


iifiz))c^ (i2))c°(fj,)S‘ 


izyiui vy2U2 




((/(^))ci(i/))ci(/x)6 


iif{z))c0{i2))c^(y.)e 


-y-yiui -izy2U2 


-/2yiui „+!^y2«2l 


g{y)dzdy. 


We subsequently note: 


((/W)c22(,2))cu(^) = ^iifiz))cn(^) - (-l)^>/(^)ci2M/r) 


= -^ifiz) - i-lY^iyfiz)iz-{-iy^yfiz)g+{-iy^+^^yi^fiz),^y) 
= \ifi^) - i-iy^i'fiz)v - {-iy^yf{z)g + {-iy^+y^vyf{z)yiy) 

((/('^))c2l (/2 ))c'? 2 (ly) /(■^)c[2l’■^ 2 ] 5 


SO the order in which we take (anti)commuting parts is irrelevant when {y, v} = 0. 
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Distributing the exponentials in UZ2U2, which are the second factors, we get four terms. Let us name 
them li with i G {1,2, 3,4}. The I3 term has the same sign as Ji in the exponential, and we 

therefore expect it to simplify with the Ii term. Indeed, we get: 

47r dR2 Jr2 

+ ((/W)cH.))cO(;.)e-™^“^e+‘'^=“^ + g{y)dzdy 

47r Jr2 Jr2 

+((/(^)),.(.)),o(^)e'^^^“^e™ + g{y)dzdy. 


If we add these two terms and combine the exponentials in y with similar signs, we notice that we have 
obtained sums {{f{z))cO(u))c°{^i) + {{f{z))cO{v))c'>-(^i) = f{z)c°{u)- In this way half of the terms cancel, and 
the following remains: 


1^1 + 1^3 = — 


47r Jr2 Jr2 

g{y)dzdy 


Here we have identified the separate qFTs and used the change of signs property m- We can subse¬ 
quently bring the root v inside the qFT of g by use of formula (HU). Let us apply the inverse transform 
to one of the constituting terms of Ii + h, for example -|- izfiz){u)T'^'''{g^^’^^){u). This gives us 

the following result: 

_p z//z/)(u)J-A‘.^(g(i4))(y)) ^ 1 (j giiA) izg^°'^'>), 

where we used Definition 12.31 of the Mustard convolution, the anticommutativity of fj, and v and formula 
(ED. We now conclude that applying the inverse transform to Ii + I3 already gives us 8 terms: 


+h) = ^ [f g + f 


-vf vg + vf *^i,u vg‘^°’^'' - vf *^,1, vg''^’'^^> + vf *^,1, vg 


( 1 . 0 ) 


di.i) 


Next, it is easily checked that the sum of I2 and I4 also simplifies: 

47r Jr2 Jr2 

+ ((/(^))ciM)cO(^)e-'^"^“^e+‘^^^“= + giy)dzdy 

47r Jr2 7r2 

+ ((/(^))cq.))c0(^)e'^"^“^e'^*'=“^ + {{/{2))cT.))cT;.)e-'^"^“^e+™] g{y)dzdy. 

The main difference is the overall minus sign of I4, which yields sums like {if{z))cO(,j))cO(fi) — {{f{z))cO(,j))c'^{fj.) 
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^(/(z))cO(^)/i. We hence obtain: 


I 2 +1 A 

47r 7 r2 Jk2 


Now we can again recognize the Mustard convolutions by applying the inverse transform. Using the 
properties (HB), (HU) and the anticommutativity of y and v we finally arrive at: 

= ^ *al,v yg - yf^°’^^ *^l,v yg^^’°'’ + 

+i2^/(04) _ ^gfiOA) ^^giOA) _ ^gfiOA) ^^^(1.0) ^ ^^^(1,1) ^ 


This results in 8 other terms. We summarize our result in the following theorem: 

Theorem 3.3. Let f and g he quaternion functions on and the left qFT with y and v anticom¬ 
muting roots of —1. The classical convolution can then be expressed as a sum of Mustard convolutions by 
the following formula: 


f*9 


= ^ / 9 + f + / *u.<' 9^^’^'’ 

- vf vg + vf *f,^y vg^'^'°'> + vf *f,^y vg^'^'^'> 

- 9-9 - 99^°’^^ + 9f^°’^'’ *u,>' 99^^'°^ + 99^^'^^ 

+:2/r/(°4) g^yg _ ^g,f(°9) g,jyg(°9) _ ,yg,f(°9) g,vg^^^°) + vyf‘^°'^'> *^^y yirg^^'^'^ 


Remark 3.4. Note that in Theorem \3.3\ and in subsequent theorems, all terms, such as vf *^ y vgd^9)^ 
should be understood with parentheses {vf) *^^y {vg^'^'^^). They are omitted to avoid notational overload. 


4 The general case {/i, v} = a 

In the deduction of Theorem 13.31 we assumed {y, v} = a = Q. We can now repeat the steps above while 
keeping a £ K non-zero. We need modifications in several steps. First of all the (anti)commuting part of 
the exponential gets extra terms. Let us recall formula (p]) : 

^7iT = - (-l)V(cos(z 2 n 2 ) - v sin {Z2U2))y) 

= i(e—+ a(-l)>sin(z2U2)) 

= 1 (_l)tg>^^2«2 “(_l)JgM^2«2 _ “(_l)Jg-/^^2U2)_ 

This would give terms with two exponentials in yz 2 U 2 in equation (1151) . Such terms are undesirable if we 
want to express everything in Fourier transforms T^''^, but fortunately we can rewrite this expression in 
‘good’ terms (which do not depend on a) and ‘extra’ terms which can be treated separately. This is done 
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as follows: 


-fiyiui -i'Z2U2 _ ^ — i'Z2U2 ^-yiyiui J_ ^ — i^Z2U2 ^y.yiui 

-^cO(p.) e +e„i,,., e 


_ _j_ gl^Z2U2 


1 

2 ^ 

1 

2 * 


c 1 (m) 
a 


2 ® 


2 ® 


-IJ,Z2U2\ -IJ-yiUi 


+ _(g-<2Z2t*2 _ g;^Z2U2 _ _g 


® fJ,Z2U2 _|_ ®g-A122ll2^gMyi«l 


= ‘good terms’ + _[e^^^ 2 « 2 g-Myiui _ g-M 2 : 2 « 2 g-™nxi _ ^i^z 2 n 2 g-/ 2 Z 2 « 2 gTOiui] 


= ‘good terms’ + -[4sin ( 2 ; 2 W 2 ) sin (j/iMi)] 


= ‘good terms’ + a[e 


|sm(yiui)- 


a r _ 

= good terms + 7 [e 


l^Z2U2-^^^ — fiyiUi 




( 16 ) 

Also, when bringing the square roots inside the transform of g, we can no longer use formula (I12p but 
we have to use equation (USD. 

Let us first focus on the modifications needed in the calculation of the ‘good’ terms. 


4.1 Derivation of the ‘good’ terms 

For the ‘good’ terms we will show that bringing the roots inside the transform of g does not alter the 
convolution formula. The derivation of Section 13.21 can go unchanged up to the following point: 

+ + v]v)[u)T^’'^ {g^^’^^)(u) 


Here, we have to bring the roots inside the transform, which is now more complicated due to the lack of 
anticommutativity of /x and v. 

To make our future arguments clearer, we will first rearrange the summands in the equation above: 


where we have omitted the arguments in u and the factor 7r/2 to improve readability. 

Now consider for example the terms and 

Using property (113 we get for their sum: 
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We observe that the terms in a cancel out, and we obtain the same result as before. 

This is similarly the case for the other terms. Just to fix the idea let us consider for example the 
terms and 












Again the terms in a cancel. We conclude that we already have the following terms in the convolution 
formula: 

^ 9 + f *f^,u 9^^"°^ + f 

- vf *ii,u vg + vf *11,u vg^°’^'' - vf *p.,u + vf 

- *^i,u 99 - *11,1' 99^°’^^ + 9f^°’^^ *^1,11 99^^’°^ + 99^^’^^ 

vg,g - vg,g^^'^'> . 

In other words, the ‘good’ terms lead to exactly the same formula as in Theorem 13.31 except for the final 
swap of 9 and v in the last line. 


4.2 Derivation of the extra terms 


Next we have to continue the derivation of the extra terms in the equation (1161) . Fortunately, we can use 
the fact that sine is real-valued and hence commutes with quaternion-valued functions. We proceed as 
follows: 

{f * 9{x))(u) 

-—[ *g){x)dx 


7b2 

1 f 


„-ti(zi+yi)ui -ii(z2+y2)u2 


27r 


'B2 JB2 


f{z)g{y)dzdy 


=‘good terms’-f — e sin (z 2 M 2 )(/(^:))cO(iy) sin (yiui)e g[y)dzdy 

7r2 J^2 


a 

27T 


'R2 J^2 


e sin {z 2 U 2 )ifiz))cii, 2 ) sin {yiui)e'^y^'^'^g{y)dzdy 


= ‘good terms’ H- 

Stt , 


„-/iZlUl ( - VZ 2 U 2 


'b2 Js? 


'){’^f{z) + f{z)v) sin {yiUi)e g{y)dzdy 


JB2 „ .. 

= ‘good terms 
a f f 


^-tiziui(^^-,iz2U2 _ - f{z)y)sm{yiui)e'^y^'^^g{y)dzdy 


-f 

-f 

-f 


,-flZiUi ( -VZ2U2 _ VZ2U2 


IStT 7^2 7r2 

— [ [ 

IfiTT Jg2 Jg2 

— [ [ 

IStT 7r2 7r2 

— [ [ 

IStT ./b2 ./b2 


)iif{z){e-^^y^^^ - e>^y^^^)9e-‘'y^^^g{y)dzdy 


g-A*zi«i(g-i.z2«2 - e^^y^^^)ye-''y^^^vg{y)dzdy 


„-y,ZiUi ( -VZ2U2 _ VZ2U2 


)iif{z){e-^^y^^^ - e^^y^^^)9e''y^^^g{y)dzdy 


g-A.zi«i(g-i.z2«2 _g^^^2«2)y(^)(g-Myl«l - e^^y^^^)ye'^y^^^ug{y)dzdy 


= ‘good terms’ + Ii + I 2 + h + h, 
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where we have denoted the extra terms li with i = 1,2, 3,4. 

Let us for example work out the extra terms coming from Ji. We have to apply formula ([5]) to bring 
/r inside the transform of g: 


h = 


IGtt 


/ / 

Js? Js? 


- e-''y^'^^^y)g{y)dzdy. 


2 ' 


We can now read off the extra terms coming from Ji by multiplying the products out, reading the signs 
of the exponentials and applying the inverse qFT. This is done for the first term of Ji below: 


IStt Jr 


-i/Z2«2 




vf{z)e ^‘^^'^^e‘''^^'^^gg{y)dzdy) 


( 0 , 1 ) 


= -^yf gg 

Eventually a lengthy calculation leads to the following end result for the extra terms Ii + I2 + I3 + I4' 

I -/ yyg + yf 99 + f + yf *11,^ 99^°'^^ 

+ f *11,^ 9yg''^’°^ - yf - vf 

+ yyg - yf^°’'^^ yg - *^,,2 yyg^°’^^ - yf^°'^^ *^,,2 yg^°’'^^ 




. 1 ) 


yyg 


(1.1) 


2/(0 


. 1 ) 


^,i'y9 


(1.1) 


-.2 r 


/ 9 - f *M.i' 9 


( 1 . 0 ) 


di.i) 


We summarize the result of the derivation (sum of the ‘good terms’ and the extra terms) in the 
following theorem. 

Theorem 4.1. Let f and g be quaternion functions on K^, and the left qFT with y and y roots of 
— 1, where {y,y} = a. The classical convolution can then be expressed as a sum of Mustard convolutions 
by the following formula: 


f *9 = 


1 

4 L' 


/ 9 + f 9^^'^'’ 

- yf ^5 + '^f yg^^'^"' - yf yg^^'°^ + yf yg^^'^^ 


-yf^° 

+yyf^°' 


. 1 ) 


yg 


.1) 


1) 


^99 


-yf^° 

- y^f^° 


L,v yg 


( 0 . 1 ) 




. 1 ) 


1,1^ yg 


(1.0) 




-i,u yg 


( 1 . 1 ) 


*u,^yy9<-°'^^-yyf^° 


*u,’' vyg^^’^'’ + yyf''°’^ 




a 

+ - 


-/ y^9 + *u,'' y9 + f *u,‘' + '^f *u,'' yg^^'^'’ 

+ / *M.>^ -1"/ - ^f *u,'' yg^^’^^ 


+ yyg - * m .>^ yg- f 


( 0 , 1 ) 


Ko.i) 


yyg 


( 0 , 1 ) 




(0.1) 


-t,i> yg 


( 0 . 1 ) 


-f^°'^^ *u,’' yg^^’°^ + yyg^^’^'* + yf^°’^^ yg^^'^^ 

f *u,i' 9 - f *u,>' 9^°'^^ - f 9^^'°^ + f *v,’' 


In the case of a = 0, Theorem 14.11 reduces to Theorem 13.31 which is indeed easily verified. Note also 
that the derivation and the resulting expression of / * g in terms of / g is much more complicated 
than the ‘inverse’ expression in Theorem 12.51 
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Remark 4.2. The number of terms in the right-hand side of Theorem \4. 1\ may seem large. However, in 
practical applications such as image processing the formulas will often simplify. Indeed, typical masks for 
image filtering will satisfy eertain symmetric properties. If e.g. = g the mask g is an even function 

in the second variable, and the term in vanishes completely. 


5 Consequences 

A first immediate consequence is that we can now easily compute the quaternion spectrum of a classical 
convolution. In the case of anticommuting roots this is formulated as follows. 

Theorem 5.1. Let f and g be quaternion functions on and IF = F^''' the left qFT with g, and v 
anticommuting roots of—1. The spectrum of the classical convolution can then be expressed as follows: 

Hf*9) 

= I [j-(/).F( 5) +.A(/).A(g(°’i)) + 

- F{yf )F(vg) + F{yf)F(vg^^^^^) - F(yf)F{vg^^^^'^) + F(vf)F(vg^^^^^) 

- Fifif^°’^y)Fit,g) - F{^if^°'^^)Fi^ig^^’^y) + Fifif^°’^^)F{pg^^’°^) + F{^if^°’^^)Fi^^g^^'^^) 
+Fiiypif^°’^'>)F{pug) - F{upf^°’^^)Fiyiug^°’^'>) - Filypif^°’^'>)Fi^^lyg^^’°^) + Fiuyif^°'^^ Fipug^^’^'>) . 

This result should be compared with the convolution theorem for the qFT obtained in [7]. The 
advantage of our result is that it gives much better insight into the nature of the spectrum: it can be 
explained by suitable reflections of the original functions and multiplication by the roots g, and v. A 
more complete discussion of how to interpret the qFT spectrum, in the context of color image processing, 
is postponed to the follow-up paper [35]. 

In data and signal processing the convolution has its counterpart in cross-correlation. It is a measure 
of how strongly two signals are correlated. Classically it is defined as follows: 

Definition 5.2. The classical cross-correlation of two complex functions f and g defined on is given 
by the following integral: 

f-*^9{y)-=[ f*{x)g{x + y)dx. 

The complex conjugation (/* is used to distinguish this from the quaternion conjugation /) in this 
definition is needed to ensure a relation between the autocorrelation function /*/ and the power spectrum 
of a signal (as required by the Wiener-Khintchine theorem). Notice that even in the complex case the 
cross-correlation ★ is not commutative, in contrast with the convolution *. Because of this conjugation 
operation there are several non-equivalent definitions of the cross correlation possible. The conjugation 
operation is not relevant for the derivation of our theorem and can be easily added. That is why we stick 
with the following definition: 

Definition 5.3. The quaternion cross-correlation of two quaternion functions on / and g, is given 
by the following integral: 

f*9{y)-= f{x)g{x + y)dx. 

With this definition we can easily observe that we can state the correlation in terms of a convolution: 

f*9 = *9- 

This allows to obtain results for the spectrum by using the theorems derived in this paper. As an 
example we will state the correlation theorem corresponding to Theorem 13.31 
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Theorem 5.4. Let f and g he quaternion functions on and T = the left qFT with fj, and v 

anticommuting roots of —1. The spectrum of the correlation can then he expressed as a sum of products 
of the individual spectra hy the following formula: 


6 Conclusions and outlook 

In this paper we have established formulas to express the classical convolution product of two quaternion 
functions as a linear combination of Mustard convolutions for the left qFT. The computations leading 
to this result can be adapted without additional difficulty to other versions of the qFT. In a subsequent 
paper we plan to show how to apply our results to the design of quaternion filters for color images. 
We also plan to extend our results to higher dimension, where the role of the qFT is taken over by a 
hypercomplex Fourier transform defined using an arbitrary number of roots of -1. 
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